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ABSTRACT 

The relationships between supermassive black holes and the properties of their as- 
sociated dark-matter halos imply that outflows from accreting black holes provide a 
feedback mechanism regulating galaxy formation. Accreting black holes with weak 
or undetectable radio jets (radio-quiet quasars) outnumber those with powerful jets 
(radio-loud quasars) by a factor ~ 10 — 100, so powerful-jet outflows are often ne- 
glected. However, whenever powerful jets are triggered, there is a dramatic (factor 
> 100) step-function increase in the efficiency of feedback. We use a feedback model, 
together with the measured space density of flat-spectrum radio- loud quasars, to show 
that a powerful-jet episode probably occurred in every protocluster in the Universe. Be- 
fore jet triggering, there was time for gravitational collapse to create many (~ 10 — 100) 
surrounding protogalaxies massive enough to host radio-quiet quasars. After trigger- 
ing, the powerful jet pushes back and heats ionized gas so that it cannot fall onto these 
protogalaxies and cool. Once neutral/molecular gas reservoirs become exhausted, there 
is a synchronized shut down in both star-formation and black-hole activity throughout 
the protocluster. These considerations imply that radio-loud quasars have a profound 
influence on the evolution of all the galaxies seen in clusters today. 
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1 INTRODUCTION 

Powerful-radio-jet activity is regarded as a useful tracer of 
large-scale structure in the distant Universe (e.g. Miley et al. 
2004), but with the exception of a few studies (e.g. Gopal- 
Krishna & Wiita 2001; Rawlings 2003), its influence on 
galaxy formation is typically ignored. This is because, even 
in the young Universe where supermassive black holes have 
high accretion rates and are visible as quasars, only a few 
per cent of accreting black holes develop powerful radio jets 
(e.g. Goldschmidt et al. 1999). It is now accepted, however, 
that quasars do play a key role in galaxy formation because 
of the existence of remarkably tight correlations between the 
masses of black holes and the properties of their associated 
dark-matter halos such as velocity dispersion (e.g. Ferrarese 
2002). 

Such correlations are most easily understood in terms of 
'feedback models' (Silk & Rees 1998; Fabian 1999) in which 
the mechanical power emerging from radio-quiet accreting 
black holes injects energy into the gaseous component of the 
young galaxy over and above that which it has acquired by 
gravitational collapse. In Sec. 2 we adopt a feedback model 
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for radio-loud accreting black holes and in Sec. 3 we use this 
model, together with hierarchical clustering theory, to com- 
pare the predicted cosmic evolution in the comoving space 
density of radio sources with the evolution observed. We 
reach some new conclusions concerning the cosmological im- 
portance of radio sources in Sec. 4. 

The convention for spectral index a is that flux density 
S v tx v~ a , where v is the observing frequency, and the radio 
luminosity function (RLF), the comoving space density of 
sources per (base 10) logarithmic interval of 1.4-GHz radio 
luminosity L1.4, is assumed proportional to L^f. We assume 
throughout a low-density, A-dominated Universe in which 
h = tfo/(100 km s" 1 Mpc" 1 ) = 0.7; Q m = 0.3; fi A = 0.7; 
fib = 0.04, cr g h -i = 0.9, and n aca i ar = 1. 



2 FEEDBACK DUE TO POWERFUL RADIO 
JETS 

Here, we adopt a form of feedback model in which pow- 
erful (radio-loud quasar) jets deliver some fraction f c m c 
of their mechanical power Q to ionized gas which, prior 
to the powerful-jet episode, is bound to a number iVhaio 
of dark-matter-dominated halos (one hosting the radio-jet- 
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producing central engine, the others in the surrounding pro- 

tocluster), each of mass Mhaio and velocity dispersion a. 

We assume this injection of energy occurs over a timescale 

tiife, and that the energy is deposited as thermal energy in 

the gas (with negligible radiative losses; see Rawlings 2003), 

and that each halo has just virialized at redshift z, so that 
I i 

a oc M h 3 alo (1 + z)2 (Somerville & Primack 1999). We con- 
sider the critical point at which sufficient mechanical energy 
is delivered to the gas so that it just becomes unbound from 
each of the -/V ha io surrounding halos. This yields a scaling 
relation 
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where / gas is the fraction of the total (dark-matter- 
dominated) mass in gas, and k\ is a constant of order unity. 
The choices of normalizing constants for each variable will 
be explained in Sec. 3. 

It is clear from Equation 1 that powerful-jet activity 
can remove gas not only from a single host galaxy, but also 
from a large number (iVhalo ~ 100) of surrounding galaxy- 
sized dark-matter halos. In the local Universe, powerful-jet 
activity is confined to black holes triggered in rich clusters 
of galaxies (e.g. Cygnus A and 3C 295) , and although X-ray 
observations provide ample evidence of ionized gas being 
pushed back and heated by such radio sources (e.g. Smith 
et al. 2002), the cluster potential wells are sufficiently deep 
(a ~ 1000 km s _1 ) that the gas remains in a gravitationally- 
bound intracluster medium. This will not be true in similar 
systems at much earlier times because hierarchical structure 
formation (e.g. Press & Schechter 1974; Percival, Miller & 
Peacock 2000) demands that, in the young Universe, the 
clusters of galaxies seen today were gravitationally-unbound 
protoclusters, collections of protogalaxies lacking any deep 
extended potential well, and hence any intracluster medium. 



3 COMPARISON OF PREDICTED AND 
OBSERVED COSMIC EVOLUTION OF 
RADIO SOURCES 

The hierarchical clustering of dark-matter halos and the 
background cosmology (Spergel et al. 2003) are now so well 
understood that it is possible to use the feedback model of 
Sec. 2 to make a firm prediction for the number of powerful- 
jet episodes triggered per comoving volume per unit cos- 
mic time, modulo just one key uncertainty, the fraction of 
newly-created halos that give rise to powerful-jet activity. 
Powerful-jet activity in the local Universe gives vital clues to 
the critical features of the halos: such jets emerge only from 
massive elliptical galaxies, with a relatively narrow spread 



in black hole mass Mbh ~ 10 



,9±0.5 



Mq, corresponding to 



halo velocity dispersions a ~ 290 ± 60 km s" 1 (McLure et 
al. 2004). We introduce a factor /halo to allow for the fact 
that a halo within this range of a is a necessary, but not 
sufficient, condition for powerful-jet activity. 



Some choices of normalizing constants in Equation 1 
follow from the assumed close mapping between black hole 
and dark- matter properties: Q = 1.25 x 10 40 W corresponds 
to the most powerful jets observed (Rawlings & Saunders 
1991), and is the Eddington luminosity of an Mbh ~ 10 9 Mq 
black hole; and M ha i ~ 5 x 10 Mq is the mass of a dark- 
matter halo with a — 290 km s _1 , collapsing at z ~ 2.5 
(Somerville & Primack 1999). Other choices were motivated 
as follows: theory demands f c m c ~ 0.5 (e.g. Bicknell et al. 
1997); / gas = 0.13 corresponds to the ratio of baryons to 
dark matter, assuming that baryons in forms other than hot 
gas can be neglected; and tnf e = 4 x 10 7 yr is (for an as- 
sumed quasar accretion efficiency ~ 0.1) the mass-doubling 
timescale for Eddington-limited growth of a black hole, con- 
sistent with lower-limits on radio source lifetimes derived 
from the observed linear sizes of powerful radio galaxies 
(e.g. Kaiser, Dennett-Thorpe & Alexander 1997; Blundell 
& Rawlings 1999). 

Powerful jets are liable only to be triggered when two 
supermassive black holes coalesce (Wilson & Colbert 1995), 
and since this can only happen as the result of a major 
merger (a special class of 'halo creation event'), we have 
used hierarchical structure formation theory (Percival et al. 
2000) to predict how the trigger rate of powerful jets depends 
on cosmic epoch (Fig. 1). Note that episodes of powerful-jet 
activity are predicted to occur over a wide range of cosmic 
epochs, persisting at some level throughout the later stages 
of the 'epoch of reionization' during which a partially ion- 
ized Universe at z ~ 15 (Spergel et al. 2003) becomes fully 
ionized by z ~ 6 (Becker et al. 2001). 

Predicting the radio emission from powerful-jet episodes 
is, in general, an extremely complicated function of prop- 
erties like the time after the jet-triggering event and the 
gaseous environment, as well as observational choices like 
frequency. We therefore focus on episodes whose obser- 
vational manifestation will be flat-spectrum (a ~ 0) ra- 
dio emission which arises when the jets happen to be 
favourably oriented, i.e. within a 'beaming angle' covering 
a sky fraction /beam ~ 0.01 (see Jarvis & Rawlings 2000). 
Such 'Doppler-boosted' emission arises from synchrotron- 
self-absorbed knots at the base of the jet, and the ratio of 
its luminosity to Q should be fairly constant throughout the 
lifetime tuf c of the radio source (and relatively insensitive 
to environment and redshift effects) whereas, in contrast, 
contributions to L1.4 from extended structures will, owing 
to inverse-Compton cooling and other effects (e.g. Kaiser et 
al. 1997; Blundell, Rawlings & Willott 1999), typically be 
a strong function of time since the jets were triggered, en- 
vironment and redshift. We use the creation rate C290 for 
a — 290 km s _1 halos (from Fig. 1) to estimate, following 
Efstathiou & Rees (1988), the comoving space density $ of 
triggered flat-spectrum sources using 
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where is a constant of order unity (incorporating an as- 
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Figure 1. The creation rate of dark-matter halos as a function of redshift z for various halo velocity dispersions a (in units of km s _1 ): 
<t = 160, dot-dashed line; a = 230, dotted line; <r = 290, solid line; a = 350, dashed line. The points corresponding to v = 1,2,3 and 4 
for each a are marked, where v is the density threshold for collapse in units of the r.m.s. density fluctuation a p : v = 5 cr i t /o>(M), where 
<5crit ~ 1-7; mass M is related to a sharp k— space filter (used to smooth the density field of mean value po) by M = 6tt 2 pok^ 3 , with k s the 
cut-off value; and a and M are related by the spherical-top-hat-collapse model (Somerville & Primack 199) so that a oc iWjJalo (1 + z) 1 / 2 . 
We have estimated the halo creation rates, following Percival et al. (2000), by: (i) calculating, at each z, the value of M appropriate to the 
target a; (ii) using the Press-Schechter (PS; Press & Schechter 1974) formalism (and the <r g h -i -normalized power spectrum) to calculate 
g(z) = i/cxp —i/ 2 /2; (iii) estimating the run of creation rate with z as g(z) X h(z), where h(z) = (1 + z )[2.3+(o.036xz)/(l+o.203xz)] j g a 
fitting formula (kindly provided by W. Percival) for d5 cr it/di, where t is the cosmic time; (iv) normalizing the curves to the TV— body 
simulation data of Percival et al. (2000), and neglecting any variations of this rate with mass, over the small range of interest at z = 0. 
Note that the interpretation of 'halo creation rates' is far more straightforward at high v (say, v > 2), where very few halos are 'destroyed' 
by being subsumed in larger halos, than for v ~ 1 fluctuations (Percival et al. 2000). 



sumption that there is ~ 1 triggering event per halo as it 
evolves through an ~ 1-dex spread in Mh a i ), and /rlf (ex- 
plained fully in the caption to Fig. 2) ensures that $ is inte- 
grated over only the top dex of the flat-spectrum RLF. This 
is the regime (Fig. 2) in which the high-redshift space den- 
sity of flat-spectrum quasars is observationally constrained 
(Jarvis & Rawlings 2000). 

The crucial unknown quantity in Equation 2 is /halo- We 
have fixed this at the value (/halo ~ 0.01) delivered by re- 
quiring roughly equal normalizations for the predicted and 
measured values of $ in Fig. 2. With this normalization 
fixed, and within the considerable current uncertainties, the 
observational constraints on $ are in good agreement with 
the gradual high-redshift decline predicted by hierarchical 
structure formation theories. However, the more interest- 
ing result is that we require /halo ~ 0.01, implying that 



only ~ 1 in 100 triggering events (in the relevant halo ve- 
locity dispersion range) generate powerful-jet episodes. The 
feedback model of Equation 1 implies that each powerful- 
jet episode influences, in the young Universe, iVhaio ~ 100 
halos. By considering the conditions before and after one of 
these episodes, we will argue in Sec. 4 that the relationship 



AWo ~ — ^ 100, (3) 



established here using physical arguments based around 
Equations 1 & 2 and a measurement of the flat-spectrum 
quasar RLF, is telling us something important about the 
galaxy formation process. 
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Figure 2. Comparison of the predicted comoving space density <!> of flat-spectrum radio-loud quasars (solid line, calculated using 
Equation 2) with measured constraints on $ from existing surveys (shaded regions, 90 per cent confidence regions), from Jarvis & Rawlings 
(2000) and Dunlop & Peacock (1990). The value of /halo = 0-01 was set to obtain rough agreement in the relative normalizations. The 
value of <E> is integrated only over the top dex of the RLF (as in Jarvis & Rawlings 2000), requiring the introduction of a factor /rlf to 
account for the lower-Li.4 population. We adopt a scaling L1.4 oc Mg H (Lacy et al. 2001), so that the one-dex spread in black hole masses 
(McLure ct al. 2004) maps onto a two dex spread in L1.4. Adopting fi = 1.75 (Jarvis & Rawlings 2000), and going one dex further down 
the RLF, implies /rlp ~ 10 -1 ' 75 . Note that the model and data diverge significantly at low redshift (2 J$ 2). This is expected because it 
is well known that any simple Press-Schechter-based formalism fails to fully explain the dramatic drop in quasar activity at low redshifts. 
There are several reasons for this (e.g. Hachnclt & Rccs, 1993) but all linked to two key facts (i) that Press-Schcchtcr theory does not 
properly account for sub-halos that become part of larger collapsed systems, e.g. galaxies in virialized clusters at low redshift; and (ii) 
that mergers of baryonic systems like galaxies get strongly suppressed once they begin to inhabit larger collapsed systems in which the 
velocity dispersion a greatly exceeds the internal velocity dispersions of the 'sub- halo' galaxies (Carlberg 1990). 



4 CONCLUDING REMARKS 

In the young Universe, 'high-peak bias' (Kaiser 1984) will 
have established the regions of space destined, by present 
epochs, to become clusters of galaxies. Each of these proto- 
clusters will contain ~ 100 protogalaxies which, from Fig, 1, 
will form over a wide range of epochs and which, because of 
'high-peak bias', will be strongly clustered. Such protoclus- 
ters have been observed in the form of emission-line-emitting 
objects and 'Lyman-break' galaxies around high-redshift ra- 
dio galaxies (e.g. Miley et al. 2004). Central black holes will 
form in all the high- velocity-dispersion (a > 160 km s _1 , e.g. 
Ridgway et al. 2001) halos, presumably with masses set by 
the feedback processes common in radio-quiet quasars (Silk 
& Rees 1998); we see from Fig. 1 that these halos tend to 
be created at much earlier epochs than the higher-velocity- 



dispersion ones. Eventually, hierarchical processes will cause 
one pair of halos, each containing a supermassive black hole, 
to merge, creating a single a ~ 290 km s _1 halo with a single 
coalesced black hole, and triggering a powerful-jet episode. 
This will be a much more dramatic type of feedback event 
because the radio source injects enough energy into its sur- 
roundings that it gravitationally unbinds ionized gas associ- 
ated not only with the host galaxy, but more widely through- 
out the protocluster (see also Nath & Roychowdhury 2002, 
and refs. therein). There has been insufficient cosmic time 
for a larger dark-matter halo to form, so this process yields 
a reservoir of protocluster gas which is not yet gravitation- 
ally bound, and is now so hot that it cannot accrete back 
onto the protogalaxies. There will then be a synchronized, 
protocluster-wide shut down of activity, be it circumnuclear 
star-formation or black-hole accretion. 
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For a protogalaxy ~ 1 Mpc from the radio galaxy, fresh 
supplies of neutral/molecular gas from accretion and cool- 
ing will be shut off after ~ 3 x 10 7 yr (~ fO— times the 
light-travel time; e.g. Blundell & Rawlings f999), and star- 
formation and AGN activity will cease once the reserves of 
neutral gas have been exhausted (taking ;$ 10 yr if star- 
formation rates of ~ 1000 M© yr -1 use up neutral/molecular 
gas reservoirs of mass < 10 11 Mq, e.g. Greve et al. 2003). 
Accounting for these time lags, and the 'high-peak bias', 
it is not surprising that there seem to be significant over- 
densities of both intense starbursting systems (Stevens et 
al. 2003) and X-ray-selected AGN (Pentericci et al. 2002) 
around high-redshift radio galaxies. There is also mounting 
evidence that the amount of gas and dust in the host galaxy 
of a powerful radio source decreases as the source expands 
from various anti-correlations between different tracers of 
this material and source sizes (Baker et al. 2002; Willott et 
al. 2002; Jarvis et al. 2003). 

We conclude that powerful-jet activity represents a dra- 
matic (factor 100) step-function increase in the efficiency 
of feedback mechanisms believed to be an essential part 
of galaxy formation. Its influence is, however, far more 
widespread as sufficient energy is delivered to the proto- 
cluster environment that gas can no longer cool onto the 
~ 100 surrounding protogalaxies. Star-formation and black 
hole activity throughout the protocluster is shut down, but, 
accounting both for time before the jet-triggering event (see 
Fig. 1) and time lags after the event, there was ample oppor- 
tunity for ~ 10 — 100 of the protogalaxies to shine as radio- 
quiet quasars, and to build up their supermassive black holes 
and stellar bulges. In crude terms, this provides a natural 
explanation for the relative number of radio-quiet and radio- 
loud quasars at high redshift (e.g. Goldschmidt et al. 1999). 
The observed space density of flat-spectrum radio quasars 
(Fig. 2) predicts that a powerful-jet episode occurs in ~ 1 in 
100 protogalaxies, which is as expected if each protocluster 
experiences ~ 1 such event. 
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